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LONG-TERM  GOAL 

Much  data  in  the  ocean  is  Lagrangian  in  nature.  Its  full  use  in  ocean  prediction  could  advance 
significantly  the  Navy's  ability  both  to  predict  ocean  conditions  and  to  assess  the  optimal  strategies  for 
deploying  Lagrangian  observational  devices  and  their  associated  sensors.  The  development  of  a  fully 
operational,  integrated  data  assimilation  scheme  will  afford  such  a  naval  predictive  capacity  in  fixed 
ocean  regions  that  can  be  comprehensively  surveyed  by  Lagrangian  measuring  devices.  The  long-term 
objective  of  this  project  is  the  building  of  comprehensive  and  reliable  data  assimilation  (DA)  platfonn 
for  incorporating  Lagrangian  data  into  ocean  models.  It  will  form  the  basis  of  an  integrated  prediction 
scheme  for  the  ocean  that  can  feed  on  both  purely  Lagrangian  and  mixed  source  data. 

OBJECTIVES 

This  project  aims  to  develop  an  operational  technology  for  assimilating  Lagrangian  data.  This  new 
LaDA  platform  is  expected  to  be  particularly  effective  in  ocean  regions  where  coherent  structures  such 
as  ocean  eddies  dominate  the  circulation.  Focuses  are  put  on:  1)  extension  of  our  Lagrangian  data 
assimilation  (LaDA)  approach  to  develop  a  flexible  platform,  in  which  a  variety  of  moving  instruments 
that  may  not  be  viewed  as  Lagrangian  in  a  conventional  sense  can  be  integrated;  2)  design  of  observing 
systems  to  take  full  advantages  of  any  moving  instruments;  3)  fonnulation  of  automated  algorithms  for 
optimal  deployment  strategies  of  the  moving  instruments  to  maximize  the  infonnation  content  as 
observation  and  enhance  predictive  skill  of  the  DA  system;  and  4)  incorporation  of  dynamical  systems 
theory  to  enhance  predictive  skill,  in  particular  coherent  structures  and  tracer  fields  associated  with 
them. 

We  use  the  tenn  “Lagrangian  data”  in  the  broadest  sense,  including  the  positions  of  various  moving 
instruments  along  the  trajectories  as  well  as  time  sequence  of  the  tracer  fields.  Autonomous  vehicles 
(AVs)  that  glide  or  maneuver  in  the  ocean  are  new  types  of  moving  instruments  that  will  be 
incorporated  into  our  extended  LaDA  platform  as  a  part  of  the  control  variables.  AVs  have  come  to  be 
a  stable  and  reliable  mean  to  measure  water  properties.  The  platfonn  will  be  ideal  for  designing  and 
performing  the  autonomous  ocean  sampling  network,  adaptive  observations,  and  optimal  deployment 
plans  of  such  moving  instruments.  The  extension  will  also  have  a  capability  to  handle  high- 
dimensionality  and  nonlinearity  of  the  operational  ocean  models.  An  example  of  the  structure  of  the 
tracer  fields  is  an  ocean  color  obtained  by  the  satellites. 
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A  comprehensive  LaDA  platform  will  be  built  upon  balanced  integration  of  data  assimilation 
techniques  and  dynamical  systems  concepts.  Recent  developments  in  the  theory  and  use  of  dynamical 
systems  concepts  for  Lagrangian  analysis  will  add  a  new  dimension  to  the  LaDA  platfonn. 
Incorporating  the  geometrical  approach  will  offer  a  much  needed  template  for  observing  system  design 
for  the  LaDA  platform.  In  return,  accurate  estimation  and  forecast  of  Eulerian  flow  field  by 
assimilating  such  data  will  strengthen  further  the  reliability  of  the  observing  system  with  moving 
instruments. 

APPROACH 

Our  approach  is  to  use  the  LaDA  as  a  basis  for  the  development  of  an  operational  technology  that 
accommodates  assimilation  of  data  from  a  variety  of  measurements  by  any  types  of  moving 
instruments  including  drifters,  floats,  and  A  Vs.  Such  a  platform  will  be  developed  through  a  hierarchy 
of  models  and  methods.  Fundamental  issues  will  be  addressed  using  the  idealized  model  flows.  Basic 
tests  will  be  perfonned  on  intermediate  model  flows,  the  double-gyre,  shallow-water  model  with 
various  wind  forcings,  and  progressing  to  the  realistic  general  circulation  model  for  operational 
application  in  Gulf  of  Mexico.  We  will  investigate  to  what  extent  the  Lagrangian  data  can  be  used  to 
estimate  the  three-dimensional  flow  evolution.  We  will  improve  the  LaDA  method  to  better  deal  with 
the  chaotic  nature  of  the  Lagrangian  dynamics.  We  address  the  localization  and  inflation  issues  that 
will  help  simplify  estimates  of  error  covariance  and  develop  a  Lagrangian  Ensemble  Kalman  Filter 
(EnKF).  Ideas  from  dynamical  systems  theory  will  be  employed  in  every  stage  to  provide  a  solid, 
mathematical  template  and  foundation  while  casting  conceptual  simplicity  on  the  complex  platform. 
This  project  is  a  joint  project  by  the  two  co-PIs,  K.  Ide  (UCLA)  and  C.K.R.T.  Jones  (UNC-CH).  G. 
Verineres  is  a  postdoctoral  fellow  who  works  on  the  LaDA  for  Gulf  of  Mexico.  H.  Salman  is  a 
postdoctoral  fellow  who  works  on  the  applications  of  the  dynamical  systems  theory  to  the  LaDA. 

Liyan  Liu  is  a  graduate  student  who  researches  the  vertical  propagation  of  the  information  by  the 
LaDA. 

WORK  COMPLETED 

The  progress  is  being  made  to  push  the  Lagrangian  data  assimilation  towards  the  operational  use  for 
the  optimal  deployment  strategies. 

Issues  essential  to  the  LaDA  have  been  investigated  further.  In  particular,  due  to  the  very  nature  of  the 
Lagrangian  data  -  chaotic  dynamics  due  to  the  distinguished  hyperbolic  trajectories  that  govern  the 
flow  template  and  make  the  neighborhood  tracers  separate  at  an  exponential  rate  -  can  lead  to  a  filter 
divergence.  A  sophisticated  quality-control  method  is  being  developed  by  combining  the  dynamical 
systems  theory  and  the  expectation-maximization  algorithm. 

To  investigate  the  potential  of  surface  drifter  data  in  estimating  the  three-dimensional  ocean  flow,  the 
LaDA  is  applied  to  the  two-layer  point  vortex  model.  Vertical  propagation  of  the  observed  information 
is  examined  in  the  light  of  the  parameter  estimation. 

Towards  the  development  an  operational  strategy,  the  optimal  drifter  deployment  is  investigated 
further  using  EnKF-based  LaDA  and  the  flow  template  regulated  by  the  invariant  manifolds.  Currently 
the  scheme  is  applied  to  a  single-layer  shallow-water  ocean  model. 
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Figure  1.  Surface  Lagrangian  and  Eulerian  fields  of  the  Gulf  of  Mexico 
by  the  multi-layer  reduced  gravity  model 
[Colors  represent  the  sea  surface  height  (Eulerian  field).  The  black  lines  are  625 
drifter  trajectories  over  eight  days(Lagrangian  field).  The  green  and  red  dots 
represent  the  departure  and  arrival  points,  respectively.] 


A  multi-layer  reduced  gravity  model  of  the  GOM  is  implemented.  This  modeling  setup  has  been 
shown  to  be  the  simplest  representation  of  the  GOM  that  simulates  the  shedding  of  eddies.  It  is  a  first 
step  toward  a  detailed  representation  of  the  physics  of  the  GOM.  The  horizontal  discretization  uses  a 
curvilinear  grid.  It  has  a  horizontal  resolution  of  about  5  km  in  the  region  of  the  loop  current.  The 
inflow  and  outflow  conditions  are  the  cosine  shaped  current  at  the  Yucatan  channel  and  the  strait  of 
Florida  with  respective  width  of  160  km  and  150  km.  In  order  to  conserve  the  total  mass  in  the  Gulf, 
the  inflow  transport  matches  the  outflow  transport.  This  configuration  allows  us  to  reduce  the  number 
of  degrees  of  freedom  to  a  minimum  without  sacrificing  on  the  resolution  near  the  region  of  interest. 
The  structured  curvilinear  grid  of  the  horizontal  discrete  domain  was  created  using  elliptic  grid 
generation  techniques.  The  boundary  of  the  grid,  illustrated  in  Figure  1,  follows  the  coastline. 

For  the  assimilation  of  the  tracer  field,  we  have  added  a  functionality  to  compute  the  passive  racer  field 
as  represented  in  Figure  2. 
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Figure  2.  Lagrangian  tracer  fields  of  the  Gulf  of  Mexico  by  the  multi-layer  reduced  gravity  model 
[Colors  represent  the  normalized  concentration  of  the  tracer.  The  dashed  black 
lines  represent  the  sea  surface  height.] 


RESULTS 

By  combining  the  EM-algorithm,  the  tracer  control  scheme  can  detect  the  hyperbolic  effect  more 
reliably  and  may  lead  to  a  robust  LaDA  method. 

The  LaDA  method  can  successfully  estimate  the  highly-idealized  multi-layer  point  vortex  model 
within  the  Rossby  radius  of  deformation.  The  barotropic  mode  plays  a  key  role  as  expected. 


For  the  multi-layer  model  configuration  of  Gulf  of  Mexico,  the  current  is  confined  to  the  2  upper  most 
layers  with  the  ratio  of  transports  between  the  two  upper  layers  being  kept  constant.  The  upper  most 
layer  carries  80  %  of  the  total  volume  transport. 

For  the  development  of  theoretical  deployment  strategies,  the  problem  is  made  complicated  by  the  fact 
that  the  positions  of  Lagrangian  data  are  continuously  changing  with  time  by  the  underlying  flow  field. 
In  order  to  tackle  this  problem,  an  understanding  of  how  the  underlying  flow  field  governs  the  motion 
of  Lagrangian  data  is  required.  This  problem  is  most  naturally  addressed  using  recent  ideas  borrowed 
from  dynamical  systems  theory  to  extract  Lagrangian  coherent  structures.  These  structures  have  been 
studied  extensively  for  geophysical  flows  in  recent  years  and  are  understood  to  orchestrate  the 
evolution  and  notion  of  material  particles. 

Figure  3a  illustrates  the  locations  of  several  key  structures  that  were  identified  for  our  double  gyre  flow 
experiments.  Also  shown  are  the  deployment  sites  of  sets  of  drifters  in  relation  to  these  various 
structures.  Figure  3b  illustrates  the  convergence  of  the  data  assimilation  method  for  several  different 
sites.  For  reference  purposes,  a  case  with  no  data  assimilation  is  included.  The  figure  compares  four 
different  cases;  uniformly  released  drifters,  drifters  released  at  saddle  trajectories,  drifters  released 
within  coherent  vortices,  and  finally  a  combination  of  saddles  and  centers.  In  all  the  cases  considered, 
nine  drifters  were  used  for  the  data  assimilation.  The  results  clearly  show  that  the  uniform  launch 
produces  the  optimum  convergence  although  at  the  expense  of  having  to  release  the  drifters  in  regions 
remote  from  one  another.  In  stark  contrast,  targeting  Lagrangian  saddle  points  appears  to  produce 
equally  accurate  results  at  long  times  but  with  only  three  different  release  sites  for  the  nine  drifters. 

This  occurs  because  Lagrangian  saddle  points  can  disperse  drifters  over  a  wide  region  providing 
effective  sampling  of  the  entire  flow  field.  Therefore,  by  exploiting  the  internal  dynamics  of  the 
governing  flow  field  through  the  targeting  of  Lagrangian  coherent  structures,  we  are  able  to  optimize 
the  perfonnance  of  our  method  through  the  optimal  deployment  of  drifters. 

The  LaDA-EnKF  works  effectively  for  ocean  data  assimilation  using  a  shallow- water  model  (Figure  2) 
by  applying  the  localization  radius.  An  advantage  of  this  method  over  other  existing  methods  is  that 
the  assimilation  time  interval  can  be  as  long  as  the  Lagrangian  correlation  time  scale. 


Figure  3.  Results  for  directed  drifter  deployment  strategy. 

[a)  Lagrangian  coherent  structures  superimposed  on  the  height  field 
at  day  60.  B)  Error  in  the  height  field.] 


Judicious  computation  of  the  error  correlation  and  innovative  construction  method  for  error  correlation 
are  crucial  in  the  coastal  data  assimilation  system  using  the  variational  approach. 

IMPACT/APPLICATION 

Observation  of  the  chaotic  dynamical  system  can  be  assimilated  effectively  as  the  quality  control 
scheme. 

The  LaDA-EnKF  will  enable  us  to  assimilate  drifter  data  into  realistic  ocean  model. 

The  flow-template  approach  shows  a  promising  potential  to  the  effective  design  of  the  optimal  drifter 
deployment  strategies. 
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